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The excitation spectrum of the density collective oscillations is computed for multi-component 
molecular mixtures both with Coulomb and (repulsive) short-range interactions. Distinct sound- 
like excitations appear, governed by the short-range interaction, which differ from the ordinary 
hydrodynamic sound. The dielectric function and the structure factor are also calculated. The 
"two-sounds phenomenon" can be understood by means of the predictions of this model. 
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This paper is motivated by the "two-sounds anomaly" persistently reported over the years in water, either in normal 
conditions or undercooled, (author?) jl|-(author?) Q as well as in other liquid molecular mixtures. Inelastic neutron. 



X-ray, Brillouin and, more recently, ultraviolet scattering, either in ordinary or in heavy water, seem to indicate an 
^ : additional, faster, higher-frequency sound, propagating with velocity ~ 3000to/s up to intermediate wavevectors 
d ' (mean inter- molecular distance in water is ~ 3 A) , beside the ordinary hydrodynamic sound propagating with velocity 
~ 1500to/s. a dispersionless mode (~ lO^^s"^) was also reported sometimes[3|,0| (as well as no additional soundQ). 
The phenomenon is also documented both by simulations of molecular dynamics and experimental data in binary 
' mixtures with large mass difference (metallic alloys, rare-gas mixtures). |l8]-|17| 

CN 

We show herein that such a "two-sounds anomaly" may appear in interacting molecular systems with (repulsive) 
short-range interaction. Such a model could reasonably be related to liquid water (or other physical systems as those 
indicated above). The velocity of the sound-like excitations is independent of temperature, in contrast with the velocity 
• of the hydrodynamic sound which is governed by the adiabatic compressibility, and thus temperature-dependent. In 
[ addition, the plasma-Hke branch of the spectrum due to the Coulomb interaction may appear as another sound-like 
I ■ mode for shorter wavelengths and weak Coulomb coupling. We report here also the computation of the dielectric 
[ function and the structure factor within such a model. 

c/2 , We start with the well-known representation of the particle density 
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for a collection of N particles enclosed in volume V, where denotes the position of the z-th particle. We consider 
Q a small displacement — > -f u(ri) in these positions, as given by a displacement field u(ri), such that the particle 
Q [ density becomes 



^(^) = ^ E ^'^'^ E ^-^''[■-'+"('■''1 - ^ ^ e^^'^ [1 - + •••] (2) 
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• for qu(ri) <C 1. Now we employ a Fourier representation 
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as well as the well-known random-phase approximation 

^e*(q-q')ri (4) 



n{r) = n - in^ ^ e*'i""qu(q) , (5) 
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where n = N jV is the particle density. By comparing equations ^ and ([5]), we can see that the small change in the 
density can be represented as 

n(r) — n = 5n{r) = —ndivu{r) , (6) 
and its Fourier transform Sn{q) = — mqu(q). 

We apply this displacement-field approach to a multi-component molecular mixture consisting of several species 
labelled by i, each with Ni particles in volume V, mass and electric charge ezi, where — e is the electron charge 
and Zi is a reduced effective charge, interacting through Coulomb potentials Lpij and short range potentials Xij- The 
mixture is subjected to the neutrality condition X^i^-i-^i = where Ui = Ni/V is the particle density of the z-th 
species. We consider elementary excitations of the particle density, whose interaction energy is given by 

^^\T.I [^^^(^ - + ^^^(^ - '^')] 5n.,iy)5n,{v') , (7) 

where ipij — e^ZiZj/ |r — r'| and 5ni(v) denotes a small density disturbance which preserves the neutrality. According 
to equation ((6]) it can be represented as 5ni — —riidivui, where is the displacement field. We use the Fourier 
transforms 
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where N — ^^Ni is the total number of particles, (p{r) = e^/r and (p{q) = (p{q) = Aire^/q^. A similar Fourier 
transform is employed for the displacement field u^, which leads to (5ni(q) = — miqui(q). We can see that only 
the longitudinal components Ui(q) of the displacement field are relevant, so we may write Ui(q) — (q/q)wi(q), 
(5rii(q) = —iqui{q), with 5n* {—q)=Sni{q) , u*(— q) — Ui(q) and q) — —Ui{q). Making use of the Fourier 
transforms introduced above, the interaction U given by equation ^ can be written as 

yq 

where fijiq) — ZiZjip{q) and n = N/V is the total density of particles. We assume a weak g-dependence of Xijil), as 
for short-range potentials. 

Similarly, the kinetic energy associated with the coordinates Ui is given by 

T = --^y^ mjn.,Ui(q)Mi(-q) . (10) 
«q 

In addition, we introduce an external field (j){r), coupled to the electrical charges, which gives rise to the interaction 

V^~i-y^niZiq(f>{q)ui{-q). (11) 

iq 

The equations of motion corresponding to the lagrangian L = T — U ~ V are given by 

mjUj + Aire^ Zj ZjUjUj + q^ Xij''^jUj — —iqezjcj) , (12) 

3 3 

where we dropped out the argument q in Ui(q) and (/'(q) and neglect the weak q-dependence of Xijil) ~ Xij- In 
order to simplify these equations we take the same (repulsive) short-range potentials for all species, Xij = X > 0, and 
analyze first the homogeneous system of equations given by lfT2|) . We introduce the notations a = 47re^, b — g^x. 
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- ZiUiUi , y = - . (14) 
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Making use of these notations, the homogeneous system of equations (|T2l) can be written as 

(-w^ + aSi) X + bS^y = , 



In addition, we have 



(15) 

aSsx + {-uj"^ + bS2) y = ■ 



2 anzi bn 
Lu Ui — X H y . (16) 
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The spectrum of frequencies uj of the system of equations (fTSl can be obtained straightforwardly. It is given by 

1 
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Si + bS2 ± \ a^Sf + 2ab (251 - S1S2) + b^S. 



(17) 



The wi-branch in equation (fT7|) (corresponding to the plus sign) represents the plasmonic excitations. In the long 
wavelength limit it reads 



iol = aSi + bS'i/Si ^ujI + bS'i/Si ,9^0, (18) 
where ujp, given by 



c,2^a5i = ^ , (19) 
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is the plasma frequency. For shorter wavelengths the wi -branch approaches an asymptote given by 

ujl-^bS2^aSllS2 , q^oo. (20) 

The Ci;2-branch in equation (fTT]) (corresponding to the minus sign) represents sound-Hke excitations. In the long 
wavelength limit it is given by 

a;| = {S2 - Si/ Si) b = vy , g ^ , (21) 

where 



vs = ^{S2-Si/Si)x (22) 

is the corresponding sound velocity. We can see easily, by applying the Schwarz-Cauchy inequality to the vectors 
a-i = \/ni/mi and bi — Zi^fuifmi, that is always positive ((^2 — S\IS\) > 0). For shorter wavelengths the 
cj2-branch of the spectrum approaches an horizontal asymptote given by 

c^l ~ (1 - Sl/SiS2) , g ^ 00 . (23) 

In the limit of vanishing Coulomb coupling (a 0) the sound-branch of the spectrum becomes = bS2 = v^q^ , 
where 



xS2^xY.^ri,/m,) , (24) 



an expression which holds also for the same mass rrii = m for all particles (one component), due to the neutrality 
condition {S3 = 0). 

The above elementary excitations, which are governed by interaction, are non-equilibrium collective modes which 
might be termed density "kinetic" modes. [l^ The sound-like excitations (ci;2-branch in equation lfT7|)) may be called 
"densitons", in order to distinguish them from plasmons (wi-branch in equation l|17p ) and from the ordinary sound. 
They may correspond to the density collective modes suggested by Zwanzig for classical Hquids.p3| We emphasize 
that these sound-Hke excitations are distinct from the ordinary hydrodynamic sound. 
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Figure 1: Spectrum of density excitations given by equation ^ for the O"^^ - pi asma. 

Indeed, the interaction corresponding to the latter can be written as 

U^^ [dr [divuir)f ^--Ly- q^u{cDu{-q) , (25) 
■' q 

where k = —{\/V){dV/dp)s is the adiabatic compressibiUty {p denotes the pressure and S stands for entropy). The 
above equation is derived by making use of the change 6V = —V{5n/n) = Vdivu in volume. We emphasize that for 
thermodynamic equilibrium we have only one displacement field u(r). Equation l(25|) together with the kinetic energy 
given by equation (flOl) for Ui{q} = w(q) leads to the sound branch ojq = VqQ^, corresponding to the ordinary sound 
propagating with a velocity vo given by 




For rrii — m (one component) the above equation gives the well-known velocity vq — 1/ nnm of the ordinary sound. 
As it is well-known, it has a slight temperature dependence, through the compressibility, in contrast with the velocity 
Vs given above for the sound-like excitations. For an electrically neutral multi-component mixture it can be shown 
easily that v1/vq > rt^xK. 

If we apply equations ((2^ and ([26]) to both ordinary and heavy water (one component, neutral molecule), and assume 
that interaction x and, respectively, the compressibility n are the same for the two kinds of water, we can see that the 
two sound velocities Vs and vq exhibit a slight isotopic effect, while their ratio Vs/vq = ny'xK does not exhibit such 
an isotopic effect, in agreement with experimental data. In this case we may take vq = 1500to/s and Vs = 3000m/s 

from experimental data and get the interaction parameter x — SOeF • A (for a mean inter-molecular spacing 

. A similar picture, given by equations (|24l) and l|26p . may apply to rare-gas mixtures, while for metallic alloys the 

Coulomb coupling must be taken into account (and equation l[22|) employed). 

If we assume the existence of a dispersionless mode in water, then we may consider that water molecule is dissociated 
to some extent, and its components have an electric charge, such that the plasmonic mode given by equation lfT9|) 
can be identified with such a dispersionless mode. Various models of dissociation of the water molecule are known, 
like OH^ — iJ+ or OH^ — H^O^ . In all cases a certain mobility of the (hydrogen) cations and 0~ (oxygen) 
anions is implied. We assume here that the dynamics of liquid water has a plasma-like component consisting of H^^ 
cations with density 2n and mass m (proton mass) and 0~^^ anions with density n and mass M — 16m, where n 
is the density of water. The excitation spectrum given by equations lfT7|f or such an 0~^^ — plasma is shown 
in Fig. 1. Taking ui = 10^^s~^{~ 5meV) of the dispersionless mode0|,Q,[3| as the plasma frequency uip given by 
u)p = 167me^z^//i (equation (fT9|) ). where fi = 2mM/{M + 2m) is the reduced mass, we get z ~ 3 x 10~^. The 
velocity of the hydrodynamic sound is given by = 1/ \/ Kn{M + 2m) according to equation l(26|) and the velocity 
of the sound-like excitations is given by Vg — y^9nx/{M + 2m) from equation (|22l) . We can see that both velocities 



exhibit an isotopic effect, but their ratio Us/wo — Sriy'xK does not, in agreement with the experimental data. From 

Vg ~ 3000to/s we derive the interaction x — 7eT^ • A"^. Similar results are obtained for other forms of dissociation, like 
OH^ — iJ+ or OH~ — H^O'^ . In this respect, the 0~^^ — plasma model can be viewed as an average, effective 
model for various plasma components that may exist in water. 

According to equation i(20|) . for shorter wavelengths the o^i-branch approaches an asymptote given by uj^ ^ 65*2 + 
aS'l/S'2. In the limit of weak Coulomb coupling this c^i-branch may appear as an "anomalous" sound given by 

= \/hSl = Vaq , (27) 

propagating with velocity 

' (28) 



— S^l S1S2 

(which is always a positive qauntity). This additional, anomalous sound is always faster than the sound-like excitations 
propagating with velocity Vg, since 

- = , \ > 1 ■ (29) 

Vs ^1 - Si/SiS2 

It is worth noting that the molecular dynamics studies which originary predicted such a fast, anomalous sound[23| 
employed indeed a Coulomb interaction and a short-range one. We note, however, that the velocity as given by 
(|28l) does not depend on the Coulomb coupling. In the plasma model for water discussed above the ratio Va/vg is 
approximately 2 (~ y/2M/9m + 5/9), but it exhibits an isotopic effect, which does not seem to be supported by the 
experimental data. 

It is easy to derive the dielectric function in the Hmit of long wavelengths from equation lfT2|) . Indeed, for charged 
particles equation Sui = —UidivUi is equivalent with Maxwell equation divEi = AnqiSrii, where the electric field is 
given by = —AnqiriiUi and qi — ezi is the electric charge of the i-th species. It follows that the internal field is 
given by 

Eint = -47re ^ ZiHiUi (30) 

i 

We get easily this field from equations (fT2l) . 



E^nt - -^g^ 2 ' 2 (31) 



in the long wavelength limit (it is proportional to x given by equation (fT4|) ). The dielectric function is defined by 
D — eE = e(D + Eint), where D = —iqcf) is the external field (electric displacement). We get the plasma dielectric 
function 

£ = 1 - Lol/io^ , (32) 

as expected. It exhibits an absorption edge (ujp) for very low frequencies. In the static limit it is reasonable to admit 
the existence of an additional internal field of intrinsic polarizability which removes the ui — singularity. 



We pass now to the calculation of the structure factor. From equation l|T6|) we can see that the displacement Ui is 
a superposition of the two eigenvectors of the system of equations lfT5|). which oscillate with eigenfrequencies wi^2, 
respectively. It follows that these coordinates are those of linear harmonic oscillators with the potential energy 
of the form miUj'^uf/2. The statistical distribution of the coordinates Ui in the classical limit is given by dw ^ 
exp{—miU}'^uf /2T)dui, where T denotes the temperature. We get the thermal averages 

T 

(uiUj) = -Sij . (33) 

On the other hand the structure factor defined by 

S{q,uj) = ^/ drdr'dt ((5n(r, t)(5n(r', 0)) e'lt""-""')-^'^* = 



N 



J di(5n(q,t)(5n(-q,0))( 



(34) 



(we leave aside the central peak) can be written as 



Nq^ f 



(35) 



Writing 

and making use of equation l[33|) we get the structure factor 



,(l)„iwit I ,,(2) iij2t 



S{q,uj) = A^Tq^ l^nf/n^m}j 



■\s{uj - wi) + ^(5(tJ - W2) 
^1 ^2 



We can see that the relevant sound contributions read 



(36) 



(37) 
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The relaxation and damping effects can be included in the above expressions of the structure factor. As it is well- 
known, they amount to representing the ^-functions by lorentzians. 

The short-range interaction x can be generalized to an interaction matrix Xij with distinct elements for each pair of 
species. In this case, the excitation spectrum of the density oscillations may exhibit multiple branches in general, for 
a multi-component mixture. In addition, it may have special features, like a dip in the plasmonic branch, or negative 
velocity for the sound-like excitations, which may indicate either an anomalous behaviour or unphysical situations, 
depending on the mutual magnitudes of the short-range potentials Xij ■ 

Now it is worthwhile commenting upon the validity of the approach presented above. If we keep higher-order terms 
in the expansion given by equation ([2]) (i.e. for moderate values of qu^) then additional interactions appear in 
equation ([9|, which leads to finite lifetimes for the density excitations. This means that for larger wavevectors q these 
excitations are not anymore well-defined excitations, as expected. Making use of equation l(33|) we can estimate the 
mean product qui for the sound- like branch as qui ~ y^T/niiVg , where the velocity Vs is given by equation ((22| . This 
gives rather small values for qm. For instance, for water we get qu ~ 0.5 (at room temperature), which shows that the 
wavevector q may take reasonable large values providing the displacement u is sufficiently small. For the plasmonic 

branch, the condition qui ^ 1 gives a cutoff wavevector ~ ^ rriiUjp/T for large Up-, for small values of the plasma 
frequency the condition becomes qUi ^ y^T/miV^ ^ 1. 

Another source of finite lifetime for the density excitations arises from the kinetic term. Indeed, under the displacement 
Tifc Tik + Ui(rifc), where r^fc is the position of the fc-th particle in the i-th species, a mixed term 



Hint = y^TOiV,fcU»(rifc) 

ik 



(39) 



appears in the kinetic term, where v^fc = r^fc is the velocity of the ifc-particle. It is easy to get an upper bound 
for this term, by using the Schwarz-Cauchy inequality. It is given by N (^rriivf/.y^'^ (rmufy^^ or, by making use of 
((33l) . VeT per particle, where e represents the mean kinetic energy (which depends on temperature, in principle). 
This estimation can be taken as an uncertainty in energy, leading to a lifetime r ~ h/VeT and a corresponding 
meanfree path A = VgT for the sound-like excitations. For wavelengths A much longer than the meanfree path, i.e. 
for wavevectors q such as g <C 1/wsT we are in the collision-like regime (w2T <C 1), and the collisions can establish the 
thermodynamic equilibrium (hydrodynamic regime) . In this case the ordinary sound can be propagated (with velocity 
vq). For q ^ 1/vsT we are in the collisionless regime, the ordinary sound is absorbed, and the non-equifibrium sound- 
like excitations ("densitons") can be propagated (with velocity Vg). Unfortunately, it is difficult to have a refiable 
estimation of the energy e, and so of the threshold wavevector qt — 1/vsT ~ ^/eT/hvs. For e = lOmeV (and 
Vs ~ 3000r7i/s, T — SOOiT) we get qt ~ O.lA , which is in a reasonable order-of-magnitude agreement with the 
experimental data.(ll]-jd|.[ll]|.(l^4l5l| It is interesting to note that if we apply this estimation to weakly-interacting 
gases, where we may take e ~ T, we get a high value of the threshold wavevector qt ~ T/hvg, since Vg is very small (the 
short-range interaction is weak). We may say that in gases there is very unlikely to exist the sound-like excitations; 



it is only the ordinary sound that exists. On the contrary, the colhsion-hke regime is quite unhkely in ordinary soHds, 
so we have there sound-Hke excitations and to a much lesser extent ordinary sound. 

Finally, we note that the collective excitations derived above contribute to the thermodynamics of liquids. Indeed, 
the free energy can be written as 

F = Fo + Fi+i^2 = -Fo+T^ln(^l-e-'""^i/^) +T^ln(l-e-'*'^^/^) , (40) 

q q 

where Fq is the free energy associated with the particle movements and <x'i_2 are given by equation (flTl) . The evaluation 
of integrals in equation (PO)) depends on the particular magnitude of the excitation spectrum, but usually the integrals 
are rapidly convergent and their contribution to the thermodynamic properties of the liquid is small. For instance, 
the sound-like contribution is approximately given by F2 ~ —TT'^VT{T/hvs)'^ /90, which is indeed a small correction 
to F(i (the latter being governed mainly by the liquid cohesion). 

In conclusion, we have shown that in interacting molecular systems there may appear sound-Hke excitations controlled 
by short-range interactions, distinct from the ordinary hydrodynamic sound. The former are non-eequilibrium excita- 
tions, while the latter appear through equilibrium, adiabatic processes. The velocity Vs of the sound-like excitations 
is independent of temperature, while the velocity vq of the ordinary sound depends on temperature, through the 
adiabatic compressibility. In order to distinguish them we propose to call the former "kinetic" modes of particle 
density, or "densitons". In addition, in the presence of Coulomb interaction, the well-known plasmonic branch is 
present in the spectrum of the density excitations, which, for shorter wavelengths and weak Coulomb coupling may 
look Hke another, anomalous, fast sound. We have shown that the "two-sounds anomaly" reported in Hquids like 
water, rare-gas mixtures, metallic alloys, etc, and documented by molecular dynamics studies, can be understood on 
this basis. 
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